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Abstract. Calculations of the statistical equilibrium of mag- 
nesium in the solar photosphere have shown that NLTE pop- 
ulations hardly affect Mg line formation in the Sun. However, 
in metal-poor dwarfs and giants the influence of electron col- 
lisions is reduced, and the ultraviolet radiation field, enhanced 
due to reduced background line opacity, results in more pro- 
nounced NLTE effects. In the photosphere of a cool star excita- 
tion and ionization due to collisions with neutral hydrogen can 
outweigh electron collisions. Analyses based on NLTE popu- 
lations lead to significantly higher Mg abundances than those 
calculated from LTE. We calculate magnesium abundances in 
10 cool dwarfs and subgiants with metallicities from —2.29 
to 0.0. The results are based on spectra of high-resolution and 
high signal-to-noise ratio. Stellar effective temperatures are de- 
rived from Balmer line profiles, surface gravities from Hippar- 
cos parallaxes and the wings of the Mg ib triplet, and metal 
abundances and microturbulence velocities are obtained from 
LTE analyses of Fe II line profiles. For stars with metallici- 
ties between —2.0 < [Fe/H] < —1.0 abundance corrections 
A[Mg/H]NLTE-LTE ^ 0.05 — 0.11 are found. As expected 
the corrections increase with decreasing metal abundance, and 
they increase slightly with decreasing surface gravity. We also 
calculate the statistical equilibrium of magnesium for series of 
model atmospheres with different stellar parameters and find 
that A[Mg/H]NLTE-LTE increases with effective temperature 
between 5200 and 6500 K. For extremely metal-poor stars the 
abundance corrections approach A[Mg/H]NLTE-LTE ^ 0.23 
at [Fe/H] - -3.0. 
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1. Introduction 

Abundance patterns of a-elements, such as Mg, Si and Ca, in 
stars with different metallicities are important for understand- 
ing the chemical evolution of Galaxy. The first stellar genera- 
Send offprint requests to: G. Zhao 
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vatory, Xinglong, China 



tions are supposed to have produced mostly a-elements dur- 
ing massive supernova explosions of type II. Studies of a- 
element synthesis in SNe II have been undertaken by Arnett 



( |1991D , Thielemann et al. (1996|), Woosley & Weaver ( |1995[ ), 
and Nakamura et al. (1999). Due to small differences in the 



way stellar winds and semi-convection are treated and in the 
specification of the mass cut and explosion mechanism their 
predictions differ, mainly in the respective Fe yields but also in 
the pre-explosion yields of Mg. There is no question that Mg, 
in principle, is less affected by the fine-tuning of SN II explo- 
sions. Therefore, inasmuch as the products of SN II nucleosyn- 
thesis are mixed into interstellar space ^"'Mg should constitute 
a reliable reference of the early evolution of the Galaxy. Due 
to the number of strong absorption lines found in the spectra 
of even the metal-poor stars, Mg I is easier to observe than e.g. 
O I. However, it shares the disadvantage of most neutral met- 
als in the atmospheres of moderately cool stars with Mg II be- 
ing the dominant ionization stage at temperatures above 5000 
K. Consequently, neutral magnesium is sensitive to deviations 
from local thermodynamic equilibrium, particularly as its ion- 
ization balance is dominated by large photoionization cross- 
sections from the 3p ''P" state. In metal-poor stars reduced 
line-blanketing produces an increased UV radiation field that, 
in combination with the large cross-sections, leads to strong 
photoionization rates; this makes Mgl even more sensitive to 
NLTE and affects any careful abundance analysis of Mg in 
these objects. Moreover, the density of free electrons decreases 
roughly in proportion to metal abundance, and the collision 
rates become correspondingly smaller. Together this change 
of microscopic interaction processes should lead to substan- 
tially stronger deviations from LTE in metal-poor stars at op- 
tical depths logTc between —3 and 0, if the reduced electron 
collisional interaction is no t compensa ted by hydrogen colli- 
sions (Baumuller & Gehren |l996]|l997| ). 

In our present study, we use several Mg I lines chosen from 
the solar spectrum which were carefully analyzed previously 



(Zhao et al. 1998, hereafter Paper I). The standard atmospheric 



models we use in our present analysis are described in Sect. 2. 
We start with test calculations including the standard model 
atom of Paper I, and we fit the calculated NLTE equivalent 
widths for models of different stellar parameters. We determine 
Mg abundance corrections due to differences between NLTE 
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Table 1. Abundance differences between NLTE and LTE obtained by fitting calculated LTE equivalent widths of Mg I lines using 
the same parameters for various stellar model atmospheres based on our line-blanketed model grid. Results refer to A[Mg/H] = 

[Mg/H]NLTE - [Mg/H]LTE 



Toff logg [Fe/H] 457lA 4703A 4730A 5528A 571lA 7657A 8806A 8213A 11828A A[Mg/H]±a 
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3.50 
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0.120 
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3.50 
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3.50 


-3.0 


0.201 
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0.206 
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0.210 


0.212 


0.429 


0.210 


0.234±0.073 
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4.50 


0.0 


0.042 


0.041 


0.046 


0.045 


0.083 


0.100 


0.042 


0.095 


0.047 


0.060±0.025 


6500 


4.50 


-1.0 


0.077 


0.102 


0.088 


0.090 


0.123 


0.150 


0.051 


0.188 


0.072 


0.105±0.043 


6500 


4.50 


-2.0 


0.096 


0.142 


0.101 


0.120 


0.105 


0.120 


0.075 


0.178 


0.094 
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4.50 
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0.200 


0.178 
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and LTE analyses for series of model atmospheres with the dif- 
ferent stellar parameters in Sect. 3. All Mg I lines considered 
here are reproduced using standard NLTE line formation tech- 
niques with the radiative transfer solved in the Auer-Heasley 
scheme, taking into account the interaction processes between 
all levels of the Mg model atom. A previous investigation (Pa- 
per I) has demonstrated how the line profiles observed in the 
solar spectrum can be used to improve the atomic model of neu- 
tral magnesium. This model is the subject of further scrutiny in 
the light of information coming from comparison of test cal- 
culations with line profiles in the spectra of metal-poor stars. 
We then calculate magnesium abundances of 10 cool dwarfs 
and subgiants with metallicities from —2.29 to 0.0 using both 
LTE and NLTE line formation techniques in Sect. 4. Our con- 
clusions are given in last section. 



2. Model atmospheres 

The statistical equilibrium calculations are performed in hori- 
zontally homogeneous LTE model atmospheres in hydrostatic 
equilibrium. We account for metallic and molecular UV line 
absorption using Kurucz' (1992) opacity distribution functions 
(ODF) interpolated for the proper stellar metallicities. In or- 
der to use the meteoritic iron abundance Xogepe = 7.51 as the 
solar reference, we scaled Kurucz' ODF by —0.16 dex to put 
his iron opacity on the proper scale. Bound-free opacities were 



computed with solar abundances from Holweger (1979), scaled 
by the stellar metallicities, and in case of [Fe/H] < —0.6 cor- 
rected for a general enhancement of a-element abundances of 
O, Mg, and Si with [a/Fe] = 0.4. Since Mg is the most impor- 
tant electron donator among the a-elements this approximation 
follows the results found by Fuhrmann et al. (1995). Convec- 
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Fig. 1. NLTE corrections of the magnesium abundances versus 
[Fe/H] for the model grids, where A[Mg/H] = [Mg/H]NLTE - 
[Mg/H]LTE- Top: Toff = 5780 K, log g = 3.50 and 4.44. Bottom: 
log 5 = 4.50 (4.44 for T^s = 5780 K), T^g = 5200, 5780, and 
6500 K. 
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lion is taken into account parametrically with a mixing-length 
of 0.5 pressure scale heights (see Fuhrmann et al. 1993| ). 

The reason for adopting this model as a standard is that 
we can easily use the model to calculate stellar atmospheres 
differentially with full physics included. On the other hand, the 
abundance differences between the LTE and NLTE calculations 
have nearly the same value independent of the type of model 
atmospheres adopted, even though we had to adjust some line 
parameters (e.g. logCe, log g/ ) in order to obtain a better line 
fit (cf. Paper I). 

3. NLTE line formation 

3.1. The standard atomic model 

The standard atomic model we used in this analysis is basically 
the same as that used for the solar calculations {cf. Paper I). 
The neutral magnesium model includes all levels n£ up to n 
= 9 and £ < n — 1 resulting in a total of 83 Mg I terms plus 
the doublet ground state of Mg II. All e nergie s were taken from 
the compilation of Martin & Zalubas ( 1980 ) except for those 
terms with £ > 5 which we re obt ained using the polarization 
formula of Chang & Noyes ( 1983| ). Fine structure splitting has 
been neglected. This model is nearly the same as that used by 
Carlsson et al. ( 1992 ). All terms are coupled by radiative and 
collisional interactions as described in more detail in Paper I. 

We have shown in Paper I that collisional interaction with 
neutral hydrogen atoms will not affect the results obtained for 
the solar Mg I lines except those in the infrared. However, 
since neutral hydrogen collision rates do not depend on metal 
abundance, they can maintain an important influence in metal- 
poor stellar atmospheres and raise the total collision rates to 
high enough efficiency at which they can even compensate 
the strongly increased photoionization. As outlined in Paper I 
the hydrogen collisions are based on the formula of Drawin 
(1969). It is most often used in a form due to Steenbock & 



on the method of complete linearization as described by Auer 
& Heasley (1976). The calculation includes all radiative line 



Holweger (1984) although the validity of this formula cannot 



be judged. We have used it as an order of magnitude estimate 
and carefully investigated the influence of hydrogen collisions 
on our atomic model by fitting all available Mg I lines modi- 
fying Drawin's formula by a scaling factor For Mg I it turns 
out (as was the case for aluminium) that this factor should 
vary exponentially with upper level excitation energy £"„ (in 
eV), Sh = 1000 exp{— n£'n/2}. This was determined in a 
fully empirical manner recomputing the complete NLTE line 
formation with statistical equilibrium equations including the 
scaled hydrogen collision rates, and it enabled us to fit solar 
lines of different excitation energies. The decrease of hydrogen 
collision cross-sections wit h excitation energy is also in rough 
agreement with Kaulakys' ( 1985 1986 ) prediction for Rydberg 
transitions. 

3.2. Statistical equilibrium calculations for different stellar 
parameters 

The statist ical eq uilibrium is calculat ed usi ng the DETAIL code 
(Giddings 1981; Butler & Giddings 1985) in a version based 



transitions, most of them represented by Doppler profiles; 99 
lines were linearized. The Mg b lines are treated with full radia- 
tive and van der Waals damping. The linearized line transitions 
were selected from test calculations including different combi- 
nations with a preference for the stronger transitions including 
the £ = n — 1 levels. Adding more transitions did not change 
the results. The bound-free transitions of the lowest 22 levels 
were linearized, too. 

In a star such as the Sun the flux in the ultraviolet spectral 
region is determined to a large part by opacities due to metal 
line absorption. As in our analysis of the solar Mg spectr um 
we use the opacity distribution functions of Kurucz (1992) to 
represent this opacity. In these ODE data single line opacities 
in small frequency intervals are represented by superlines; con- 
sequently, the ODE opacity is not identical to that required at 
a specific position in frequency space. Eor the calculation of 
a model atmosphere this simplification is a sufficient approxi- 
mation, but for NLTE line formation the actual radiation field 
across a line transition or an ionization continuum is important 
for the determination of the statistical equilibrium of an atom. 
Eor bound-free transitions the exact position of the absorbing 
lines is less important, and the use of the ODEs will be reli- 
able, provided the intervals are small enough in the frequency 
region near the ionization edge. Eor a bound-bound transition 
with its narrow line width it can be important in which part of 
the broad synthetic ODE line it is formed. We include the ad- 
ditional ODE opacity in the UV for wavelengths between 1300 
and 3860 A to allow for a realistic behaviour of the ioniza- 
tion from the ground state and the first excited level without 
affecting most of the line transitions. Only a few Mg I lines are 
found in this region allowing us to omit additional line opaci- 
ties as most of these transitions are of minor importance for the 
statistical equilibrium. 

We use the statistical equilibrium of Mg I with hydrogen 
collisions from Drawin's formula (Drawin 1969) but scaled ex- 



ponentially with upper level excitation energy as described in 
Sect. 3.1. Abundance differences between the NLTE and LTE 
calculations with various stellar models are shown in Table |l]. 
These results are obtained requiring NLTE and LTE line forma- 
tion to fit the same equivalent width. The results confirm that 
deviations of the Mg I level populations from LTE increase with 
decreasing metal abundance. The NLTE effects in solar metal- 
licity stars are almost negligible. In Eig. [l| we demonstrate the 
variation of the NLTE corrections of magnesium abundances 
versus [Ee/H] for sequences of model grids with constant sur- 
face gravity, log 5 - 4.5, and T^s = 5200, 5780 and 6500 K, 
and with constant effective temperature. Toff = 5780 K, and 
log (7 = 3.5 and 4.44. Erom these plots it is evident that the 
NLTE effects are systematically stronger for the hotter models, 
which is in accordance with the statistical equilibrium of alu- 



minium and sodium (Baumiiller & Gehren 1996; Baumiiller 
et al. 1998). The dependence upon temperature is, however. 



strongly tied to the change in ionization; as soon as Tcs drops 
below 5300 K, Mg II starts to recombine, and Mg I becomes 
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Fig. 2. Top: Magnesium abundances of 10 cool dwarfs and 
subgiants from NLTE and LTE line formation. The symbols 
of triangle and filled circle represent the LTE and NLTE re- 
sults, respectively. NLTE corrections A[Mg/Fe] of the 10 pro- 
gram stars versus [Fe/H], where A[Mg/Fe] = [Mg/Fe]NLTE - 
[Mg/Fe]LTE as a function of metal abundance (middle) and sur- 
face gravity (bottom) 



more dominant. For stars at the turnoff the temperature de- 
pendence is negligible. As expected, the strongest departures 
from LTE are found for models with high temperature and low 
metallicity. The reduction of surface gravity results in a de- 
creased efficiency of collisions by electron and hydrogen atoms 
which again leads to slightly stronger NLTE effects. 



4. Magnesium abundances in cool stars 



Our present study is based on spectroscopic data obtained by 
K. Fuhrmann and M. Pfeiffer using the fiber optics cassegrain 
echelle spectrograph (FOCES) fed by the 2.2m telescope at the 
Calar Alto Observatory during three observing runs in Septem- 
ber 1995, October 1996 and May 1997. The 1995 spectra 
were exposed to a 1024 x 1024 24^ CCD, and the resolv- 
ing power was 35,000. In October 1996 and May 1997 a 
2048 X 2048 l5^J. CCD was employed at R - 60,000. Almost 
all the stars were observed at least t wice w ith the signal-to- 
noise ratio of about 200 (Fuhrmann 1998| ). Details with re- 
spe ct to the FOCES spectrograph can be found in Pfeiffer et 
al. ( |1998| ). For HD 19445, HD 95128,HD 103095, HD 194598, 
and HD 201891 additional spectra were obtained by G. Zhao 
and H. Zhang with the Coude Echelle Spectrograph attached to 
the 2.16 meter telescope at Beijing Astronomical Observatory 
(Xinglong, China) in October 1996 and March 1997. The de- 
tector is a Tek CCD (1024x 1024 pixels of 24 /im^ each). The 
Coude Echelle Spectrograph has two beams: a blue path and a 
red path. For the blue path, the 79 grooves/mm echelle grating 
is used along with the prism as cross-disperser and a 0.5 mm 
slit, leading to a resolving power of the order of 44,000. For 
the red path, the 31.6 grooves/mm echelle grating is used along 
with the prism as cross-disperser and a 0.5 mm slit, leading to 
a resolving power of the order of 37,000. The more detailed 
technique and performance descriptions of the Coude Echelle 
Spectrograph can be found in the paper of Jiang (1996). 

The effective temperatures of our program stars were deter- 
mined from the wings of the B aimer Unes Hq, and profile 
fitting. For metal-po or star s, and Hg are also taken into ac- 
count (cf. Fuhrmann 1998 ). Fortunately, the dependence of the 
results on effective temperature is not important for most of the 
stars; only HD 6582 and HD 103095 are cool enough to be af- 
fected (see Fig. |l]). Surface gravities have been derived from 
Hipparcos data using the method given by Nissen & Schuster 
(1997). The final surface gravities are determined by compari- 
son to the strong line wings of the Mg lb triplet which gives a 
slightly higher value than that from the determination of Hip- 
parcos data. The [Fe/H] values were adopted from the deter- 
mination of Fe II lines which are insensitive to NLTE effects 



(Thevenin & Idiart 1999). The mircoturbulence velocity was 
estimated using an extended set of Fe lines. 

We note that only a few Mgl lines show a suitable line 
profile for spectrum synthesis in our spectra of sample stars. 
Therefore, four fines, namely AA 4571, 4703, 5528 and 571 1 A, 
which are often used for magnesium abundance determinations 
are selected in our present analysis. Synthetic spectra using 
either LTE or NLTE level populations were calculated using 
the line data given in Paper L The final results are presented in 
Table ^ where both the NLTE and the LTE line profile fit abun- 
dances are listed for comparison. We present the magnesium 
abundances of 10 cool dwarfs and subgiants using NLTE and 
LTE line formation in Fig. ^ the plots of the NLTE correction 
A[Mg/Fe] versus [Fe/H] and log 5 are shown, too. One of 
the most striking results is the difference between NLTE and 
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Table 2. Stellar parameters and magnesium abundances [Mg/Fe] calculated with LTE and NLTE line formation 
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Object V Tcfi logs [Fe/H] 4571 A 4703 A 5528 A 5711 A 

HD [mag] [K] [km/s] LTE NLTE LTE NLTE LTE NLTE LTE NLTE 



6582 


5.16 


5387 


4.51 


-0.83 


0.89 






0.33 


0.39 


0.36 


0.43 


0.37 


0.43 


19445 


8.06 


6016 


4.43 


-1.95 


1.35 


0.38 


0.48 


0.13 


0.36 


0.26 


0.36 






30743 


6.26 


6298 


4.09 


-0.45 


1.64 


-0.02 


0.04 






0.07 


0.14 


-0.02 


0.10 


61421 


0.37 


6470 


4.02 


-0.03 


1.95 


-0.19 


-0.14 






-0.01 


0.03 


-0.13 


-0.01 


95128 


5.05 


5892 


4.27 


-0.00 


1.01 










0.02 


0.04 


0.05 


0.12 


103095 


6.45 


5110 


4.67 


-1.35 


0.85 


0.21 


0.27 


0.26 


0.28 


0.31 


0.32 


0.25 


0.29 


140283 


7.21 


5810 


3.68 


-2.29 


1.49 










0.14 


0.27 






165401 


6.81 


5811 


4.43 


-0.39 


1.10 


0.17 


0.23 


0.17 


0.23 


0.23 


0.27 


0.18 


0.26 


194598 


8.34 


6058 


4.33 


-1.12 


1.45 


0.25 


0.33 


0.20 


0.32 


0.24 


0.32 


0.19 


0.31 


201891 


7.37 


5943 


4.28 


-1.05 


1.18 






0.31 


0.43 


0.38 


0.45 


0.29 


0.43 




Fig. 3. Synthetic flux profiles of selected Mgl lines (dashed) compared with the observed spectra of program stars (continuous). 
The synthetic line profiles refer to NLTE with hydrogen collisions scaled exponentially with excitation energy, (top left): inter- 
combination fine 3s ^S - 3p^P° at 4571 A. (top right): excited fine 3p^P° - 5d^D transition at 4703 A. (bottom left): excited 
line 3p^P° -4d^D transition at 5528 A. (bottom right): excited line 3p^P° -5siS transition at 5711 A. 



LTE abundances in the very metal-poor subgiant HD 140283 
and subdwarf HD 19445, corresponding to factors of 1.35 and 
1.39, respectively. It is again evident how NLTE abundance 
corrections increase with decreasing stellar metallicity. There 
is one star HD 103095 (Groombridge 1830) that shows a 
different behaviour in Fig. |^ (middle); due to its low effective 
temperature (^^ 5110 K, the lowest one in our sample stars) its 
ionization equilibrium is more affected by temperature than by 
NLTE effects. This is in accordance with our analysis for the 



grids of the model atmospheres in Sect. 3. The effective tem- 
perature, surface gravity, metallicity and microturbulence have 
typical eiTor of ATeff ~ 80 K, Alogg ~ 0.1 dex, A[Fe/H] - 
0.07 dex and A^t ^ 0.2 km/s. Consequently, the error estimate 
for the magnesium over hydrogen abundances is about 0.07 
dex. These error estimations are nearly the smae for LTE and 
NLTE abundance determination. Note that in most cases the 
abundance differences are well above the observational errors. 
On the other hand, the NLTE abundances based on various line 
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profile fits generally produce a satisfactory smaller standard 
deviation than that of LTE abundances. This result is perhaps 
the most important because it shows that with high resolution 
{R > 40000), high signal-to-noise spectra it is possible to 
reproduce line profiles with very high accuracy. A few profile 
fits for the different magnesium lines are given in Fig. |] to 
show the quality of the line synthesis. They are representative 
for the average fit quality obtained for the R ^ 40000 and 
60000 spectra. 



5. Conclusions 

The variation of [Mg/Fe] with the stellar mean metallicity 
[Fe/H] contains information about the chemical evolution of 
the Milky Way. Previous investigations of magnesium ba sed on 
LTE analyses by a number of researchers (e .g. Ma gain 1987 ; 



Hartmann & Gehren 1988 



1995 Fuhrmann 



Fuhrmann et al. 
1998| ) show that magnesium has a solar abundance ratio in 
disk stars ([Fe/H] > —0.5), and then [Mg/Fe] becomes con- 
stant at +0.3 ...0.4 in metal-poor stars ([Fe/H] < -0.5). As 
shown by Fuhrmann ( 1998| ), this overabundance is a distinct 
result of population membership. Our new results based on the 
NLTE analysis of 10 cool unevolved stars demonstrate that the 
magnesium overabundance with respect to iron becomes even 
larger by roughly 0. 1 dex for stars with reduced stellar metal 
abundance [Fe/H] due to photoionization dominating over col- 
lisions. The reduction of the surface gravity results in a de- 
creased efficiency of collisions by electron and hydrogen atoms 
which also leads to small NLTE effects. The strongest depar- 
tures from LTE are found for metal-poor stars with high tem- 
perature and low gravity, i.e. turnoffstais. 

Fortunately, the trend obtained from LT E abu ndance analy- 
ses of magnesium in cool stars (Fuhrmann 1998 ) is confirmed. 
The significant [Mg/Fe] ratio in metal-poor stars is found from 
NLTE spectrum analyses, with the increased level indicating 
a shift to even higher stellar masses in the parent population 
(III). Such a trend was also found by Fuhrmann et al. (1995) or 
Audouze & Silk (1995) who investigated the initial mass func- 
tion based on the a-element enrichment of population II. Re- 
cent investigations of the Fe Il/Fe I equilibrium (see Thevenin 
& Idiart 1999 ), however, seem to demonstrate that [Mg/Fe] 
may be more affected by deviations from LTE when using Fe I 
linese to derive iron abundances. We recall that in our analysis 
we used only Fe II linese which are insensitive to NLTE effects. 

In order to confirm the present results for magnesium, and 
in particular the mean abundance differences between LTE and 
NLTE it is necessary to reinvestigate a larger sample of stars 
with various metallicities based on the NLTE analysis, espe- 
cially for the metal-poor stars of the thick disk and halo. Such 
an investigation should, however, include an Fe I NLTE anal- 
ysis. Then it will be an important step improving our present 
understanding of stellar nucleosynthesis and the chemical evo- 
lution of the Galaxy. 
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